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Abstract. We have performed a large ground-based search for transiting Hot 
Jupiter planets in the outer regions of the globular clusters 47 Tucanae and u> 
Centauri. The aim was to help understand the role that environmental effects 
play on Hot Jupiter formation and survivability in globular clusters. Using the 
ANU lm telescope and a 52' x 52' field, a total of 54,000 solar-type stars were 
searched for transits in both clusters with fully tested transit-finding algorithms. 
Detailed Monte Carlo simulations were performed to model the datasets and 
calculate the expected planet yields. Seven planets were expected in 47 Tuc, 
and five in lj Cen. Despite a detailed search, no planet-like candidates were 
identified in either cluster. Combined with previous theoretical studies of planet 
survivability, and the HST null result in the core of 47 Tuc, the lack of detections 
in the uncrowded outer regions of both clusters indicates that stellar metallicity is 
the dominant factor inhibiting Hot Jupiter formation in the cluster environment. 



1. Project and Dataset Overview 



Globular clusters are excellent laboratories in which to perform a dedicated 
search for planetary transits. In a single field of view the brightest globulars 
present tens of thousands of solar- like main sequence stars, which can be simul- 
taneously analyzed for the photometric signature of an orbiting planet. These 
clusters allow study into how environment affects Hot Jupiter planet formation 
and survivability, by placing strong observational constraints on the frequency 
of Hot Jupiter planets in these non-local environments. Such constraints on 
short period planet frequencies at low metallicities are not well known, due to 
the metallicity bias of the nearby bright stars studied by ongoing radial velocity 
searches. In the case of globular clusters, are the frequencies different, and if so, 
is stellar density or metallicity more important in determining the frequencies? 

To address these issues, we have sampled the outer halos of the two bright- 
est globular clusters in the sky, 47 Tuc and u> Cen. This work comprises the 
first dedicated transit survey in u> Cen, and the first wide-field ground-based 
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sur vey of 47 Tuc, follow ing on from the HST null result in the cluster core 
by iGilliland et all (l2000Tl. F or 47 Tuc, our results have previously been pub- 
lished (jWeldrake et al.ll2005h . placing a strong constraint on the frequency of 
Hot Jupiter planets (R~1.3Rj up ) out to an orbital period of 16d. 

For co Cen, our observations placed constraints on planets with larger radius 
limits (R=1.3— 1.6Rj up ) with an upper orbital period limit of 7d. Major side- 
results incl ude the detection of va riable stars, with 100 found in 47 Tuc (69 new 
discoveries; I Weldrake et al. 2004! ) and 187 found in uj Cen (81 new discoveries; 
IWeldrake et al.ll2006al ) .In total both surveys present V and I photometry (with 
astrometry) for a quarter of a million stars. 

Our datasets were produced with the ANU lm telescope located at Siding 
Spring Observatory and the 52' x 52' field of view of the ANU Wide Field Imager. 
They are most sensitive to the cluster outer halos where photometric accuracy 
is higher due to reduced crowding. These parts of the cluster have distinct 
advantages over the cores; they have small stellar densities and are dominated 
by low mass stars. We have 56 nights of data covering both clusters, com- 
prising 2100 images, and have photometry for ^53,000 cluster main sequence 
stars (17.2<V<19.5) with which to search for transits. Our total light curve 
database was produced via differential photometry and comprises 220,000 stars 
(14.0<V<22.5). Our photometric accuracy is 2% at V=18 2, an d 4% at V=19.0 
across the field, after an application of the lTamuz et al. ( 20051 ) systematics re- 
moval package. The whole field except the inner 6' of the core was sampled 
in both clusters. Considering the detection statistics, such a large photomet- 
ric database presents an excellent opportunity to detect Hot Jupiter planets and 
statistically study their relative frequencies compared to the solar neighborhood. 



2. Transit Expectations 

A significant amount of statistical analysis was carried out on both datasets in- 
dependently, to determine the transit recoverability and expected planet yields. 
The first step involves finding the expected depth and duration parameters of 
Hot Jupiter transits as a function of V magnitude on the main sequence of each 
cluster. As the cluster stars can all be assum ed to lie at the same distance, this 
was performed by producing Yi et al. ( 20031 ) theoretical isochrones to simulate 



the stellar population of each cluster. The stellar radius was found as a function 
of V magnitude along the cluster main sequence and the resulting transit depth 
and duration determined for a hypothetical Hot Jupiter planet orbiting with an 
assumed period of 3.3d (the most typical period of planets in the solar neigh- 
borhood) and radii of 1.3Rj up and 1.6Rj up . The results can be seen in Fig. Q] 
for both 47 Tuc (left) and to Cen (right). 

The bright magnitude limit for the transit search lies at the location of 
the cluster main sequence turnoff, as stars brighter than this correspond to 
foreground stars and cluster stars of rapidly increasing radius, with increasingly 
undetectable transit depths. For our clusters this limit is V=17.2 for 47 Tuc and 
V=17.5 for the slightly more distant uj Cen. Both clusters display comparable 
photometric results, and the faint limit to the search is defined as the magnitude 
at which the photometric uncertainty becomes larger than the expected transit 
depths. Transit recoverability becomes increasingly truncated after this point. 
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Figure 1.: The expected depths (magnitudes) and durations (hours) of plane- 
tary transits as a function of V magnitude on the main sequence o f 47 Tuc 
(left) and u Cen (right) using stellar radius parameters found from lYi et all 
(2003) isochrones. The duration values assume a fixed orbital period of 3.3d, 
the most typical period of Hot Jupiter planets in the solar neighborhood. The 
dotted boxes limit the search range, with the bright limit at the location of 1% 
photometry and the lower limit where our uncertainty reaches 4%. With these 
parameters, simulated transits were produced to determine the statistical transit 
recoverability in each dataset. 

This occurs at V=19.5 for both clusters. These limits allow us to sample at 
least the upper 2 magnitudes of the cluster main sequence for transit signatures. 
Appropriate transit models can be made to match the radii of this subset of 
the cluster population, testing the transit detection algorithms and allowing a 
determination of the expected planet yields. 

3. Monte Carlo Simulations and Transit Recoverability 

A full statistical analysis was performed via extensive Monte Carlo simulations 
to determine the transit recoverability in our two datasets. Simulated transits 
were added to actual dataset light curves using the expected transit depth and 
duration values for each cluster. By using real lightcurves, our simulated transits 
contain the same photometric information as our real data; this allows us to have 
a better idea of the effects of systematics on transit detectability, and the effect 
of the dataset window function on transit visibility. 

For our Monte Carlo simulations, we consider that bins of differing photo- 
metric uncertainty correspond to stars within a certain magnitude range, and 
study each in turn. First, we consider that all stars with photometric precision 
<1% correspond to brighter stars lying purely in the foreground of the cluster, 
as cluster stars of this brightness correspond to giants. These foreground stars 
have sufficiently good photometry to permit transit detection of orbiting com- 
panions down to a radius of 0.7Rj up . The upper limit was arbitrarily chosen to 
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Figure 2.: The weighted mean statistical transit recoverability for the u Cen 
dataset. This takes into account the individual recoverabilities for differing 
planetary radii (over-plotted) as well as the number of stars as a function of 
magnitude along the main sequence. Our data displays good recoverability to 
an orbital period of ~3d for the larger radius companions. 



be 1.6Rj up . Hence, for these brightest stars, simulated transits were added to 
dataset light curves corresponding to orbiting companions of radius 0.7, 1.0, 1.3 
and 1.6Rj up , with duration corresponding to the period range of 1— 16d (for 47 
Tuc) and 1— 7d for to Cen. 

The cluster was simulated in two photometric precision bins, one with 1—2% 
photometry corresponding to the cluster upper main sequence, and one with 
2—4% corresponding to the lower main sequence, extending down to our photo- 
metric limit. As the photometry becomes worse, smaller planets are no longer 
detectable, hence for the upper main sequence bin, transits corresponding to 
orbiting companions of 1.0, 1.3 and 1.6Rj up were produced, and only 1.3Rj up 
and 1.6Rj up companions for the lower main sequence. To determine sufficiently 
the transit recoverability, over a million simulated transits were produced with 
appropriate parameters for the transit detection algorithm s to recover. 

We used two detection algorithms, the BLS metho d of lKovacs et al.l ( 20021 ) 
and the MFA algorithm of IWeldrake Sz Sackettl ( 20051 ). Both algorithms pro- 
duced very similar recoverability rates, and the MFA algorithm was used as a 
first-pass test to eliminate many false detections before BLS was used to find 
accurate period information for the resulting candidates. Any candidates which 
matched the true period of the input transit, or corresponded to an integer alias 
were classed as recovered. The recoverability of each photometry bin was deter- 
mined, and then the weighted mean transit recoverability was found taking into 
account the relative numbers of stars that lie in each of the bins. 

The resulting transit recoverabilities for the case of u Cen can be seen 
in Fig. El The recovera bility for 47 Tuc has previously been published in 



Weldrake &: Sackett] ( 20051 ). The recoverability for u Cen is good in this case out 



to an orbital period of ~3d, after which it becomes increasingly truncated due 
to the temporal properties of our dataset. The larger radius companions have 



Transits in Globular Clusters 



5 



higher recoverability, as expected due to their increased transit signal. With this 
information it is possible to determine an upper limit to the expected number 
of planets we should detect in our two datasets. 



4. Expected Number of Planets 

For any transit search, accurate calculations must be made to determine the 
expected planetary yield, in order to place meaningful statistics on the relative 
frequency of planets in the target field. For stellar clusters, this task is simplified 
somewhat, as all the target stars can be assumed to lie at the same distance. 
Theoretical isochrones can give useful information on the physical properties of 
the stars inside the search regime. Several factors must be considered to find 
the planet yield, namely the total number of stars, the occurrence frequency 
of detectable planets, the transit probabilities incorporating the radii of the 
target stars, and the weighted mean transit recoverabilities - which in our case 
incorporates residual dataset systematics and the dataset window function. 

We determine separately our sensitivity to 1.3Rj up and 1.6Rj up planets and 
assume their orbital period distribution to be the same in the clusters as it is in 
the solar neighborhood. In total, we have data for 22,000 main sequence stars for 
47 Tuc and 31,000 for uj Cen. Our search encompasses the orbital period regimes 
of both Hot Jupiter and Ve ry Hot Jupiter plan ets, and thus each type of planet 
is theoretically detectable. I Gould et al.l ( 20061 ) present a statistical analysis of 



the occurrence frequencies for both of these types of planet, calculated from 
results of radial velocity surveys and transit surveys combined. They conclude 
that Hot Jupiters and Very Hot Jupiters have occurrence frequencies of 1/310 
and 1/690 respectively. These numbers can be used to find the total possible 
number of such planets that can be associated with the stars in our sample, 
assuming these frequencies hold for our clusters. 

However, not all of these planets will transit, and the transit probability was 
determined taking into account the radii of the target stars, and the weighted 
mean radius for the stars in 47 Tuc is O.94R0, with O.96R0 for ui Cen. The 
transit probability was then found for orbital periods from 1— 16d for 47 Tuc 
and 1— 7d for uj Cen with this stellar radius in bins of 0.5d width. 

By incorporating the previously determined transit recoverability for these 
same bins, the total number of detectable planets was found by multiplying 
the factors above and summing over all bins. For 47 Tuc, we expect to find 7 
planets (3.3(7 null significance) of radius 1.3Rj up with a period range of 1— 16d. 
For uj Cen, we expect 5.3 planets (2.8a null significance) in the radius range 
1.3Rj up — 1.6Rj up with period 1— 7d. These are firm upper limits to the expected 
number. When taken together, these numbers provide a high significance result 
if no planets are found. A total of 0.7 planets are expected among the foreground 
Galactic disk stars. 



5. Search Results and Conclusions 

A vigorous search was made of all 53,000 stars with both the BLS and MFA 
transit detection algorithms. Despite a detailed examination of all candidate 
detections, we found no planetary transits around any sampled cluster star. 
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This null result has a significance of 3.3a from Pois son statistics for 47 T uc and 
2.8o" for uj Cen, and agrees with the null result of ( Gilliland et al.|[2000l ) in the 
core of 47 Tuc, indicating that the occurrence frequency of Hot Jupiter planets 
is reduced in 47 Tuc and u> Cen. 

Our expected planet numbers assume that stellar metallicity does not affect 
short period planetary frequencies, contrary to the apparent trend in the solar 
neighbourhood. By accounting for the low metallicities of the clusters using 
this trend, the expected number of planets we should find reduces essentially 
to zero. The dynamics of the clusters are not sufficiently violent to disrupt 
planetary systems, with a characteris tic lifetime of close- in planets estimated 
to be greater than the Hubble time ( Fregeau et al.l [20061 ) . Considering this, 
our null result suggests that stellar metallicity, not dynamics, is the dominant 
effect limiting the frequency of short period massive planets in globular clusters, 
and places an observational constraint on planetary frequency at such a low 
metallicity. Perhaps the low metallicity does not affect planet formation, but 
does affect planetary migration. If true, then long period planets should still 
exist in these clusters, undetectable in our work. 



6. Lupus Control Field 

A control dataset in the Lupus Galactic Plane field was observed for 53 nights in 
2005 and 2006 in order to test the observing strategies and analysis methods we 
performed on our globular cluster datasets, and compare the results. Analysis 
of ~26,000 stars with photometry better than 2% is ongoing, and has led to the 
identification of transit can didates, one of which h as excellent prospects for being 
a new Hot Jupiter planet (| Weldrake et al. 20071 ). As no such candidates were 



seen in either of our clusters, with the same methods employed, these results 
strengthen the conclusion that short period planets are far rarer in globular 
clusters than in the general Galactic field. 
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